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1 INTRODUCTION 

The purpose of this paper is to present a calculation of the string suscepti- 
bility exponent for a statistical system of dynamically triangulated random 
surfaces coupled to four species of Ising spins. The method employed con- 
sists in simulating on a computer the so-called grand-canonical ensemble of 
surfaces with the topology of a sphere. The motivation for this study is to 
contribute to a better understanding of the physics of random surfaces be- 
yond the "barrier" at c = 1 (c denotes the central charge of matter fields 
interacting with the geometry). 

As is well known, the spectacular progress in our understanding of two- 
dimensional quantum gravity, achieved during the last years, concerns models 
where c < 1. The string susceptibility exponent 7 is then given, for spherical 
topology, by the formula 

l = ^[c-l-4i^^W^^)l (1) 

derived first, in the continuum framework, by Knizhnik, Polyakov and Zamo- 
lodchikov []T|. The right-hand side of (|^) becomes complex for c > 1, which 
is clearly unphysical and indicates that something rather dramatic happens 
when one attempts to move c beyond unity. It is generally believed that this 
"decease" reflects the tachyonic instability of string models. A suggestive 
heuristic picture imputes the break-down of the Euclidean theory for c > 1 
to a condensation of gravitational singularities . 

Eq. (|lD also holds in solvable models of discrete quantum gravity, when 
the "dynamical triangulation" recipe is adopted to discretize the two- 
dimensional manifold. Nevertheless, the discrete models of Euclidean quan- 
tum gravity are well defined for any value of c. Triangulated surfaces collapse 
or degenerate into branched polymers when c — > cxd 1^,0. The indications 
of the phenomenon have been observed numerically, for c > 10, many years 
ago 1^, ^. Also, there exist theorems establishing the absence of a non-trivial 
continuum limit for classes of discrete models (for a lucid discusion of these 
issues we refer the reader to an excellent review by Ambj0rn [0). However, 
in numerical simulations, one does not observe anything special to happen as 
one crosses the famous "barrier". Moreover, an analysis of a series expansion 
of the partition function, rewritten using the matrix model formalism, has 
been carried out by Brezin and Hikami P]. They also did not find any sign 



of pathology at c > 1. This situation has motivated a series of numerical 
simulations, with the aim to achieve a better understanding of the quantum 
gravity coupled to matter with central charge larger than, but close to unity. 
As these studies were restricted to the microcanonical ensemble ||^, [1^, |ll], |1^ 



we have decided to use the grand-canonical ensemble to determine directly 
the particularly interesting susceptibility exponent, for the model with four 
species of Ising spins (c=2). In parallel, a novel method of measuring this 



exponent has been developed in Copenhagen [13]. When completing this 



work, we received the very recent ref. [|J], where 7 is also measured for the 



model we have been working with. Our calculation constitutes an indepen- 
dent check of these results, obtained with a more traditional method. In 
the last section of this paper we also conjecture a specific scenario for c > 1 
models. 



2 THE NUMERICAL EXPERIMENT 

2.1 The method 

The model is defined by the partition function 

Z(^,/3) = ^iy(T)e-^^(5:exp[/3^a.a,]] , (2) 

T,N \[a] <ij> J 

where T refers to a specific triangulation and W{T) is the corresponding 
symmetry factor. Further, A^ is the number of vertices and n is the number 
of species of Ising spins a. The cosmological constant g and the spin cou- 
pling f3 are the only parameters of the model. Only nearest neighbour spins 
are coupled. As in the early ref. [J15l the spins live on the vertices of the 



triangulated surface (and not the dual one). 
The partition function can be rewritten as 



^e-^^Z(Ar) (3) 



N 



For large A^, one expects 



N^Z{N) ~ A^^ e^^"-^ , (4) 



where 7 is the susceptibihty exponent and Qcr is the (non-universal) critical 
cosniological constant. 

A single step of the grand-canonical algorithm [^, ^ changes A^ by ±1. 
The code employed in this work to update the geometry is essentially identical 
to the one used and described in detail in ref. [^. The only modification 
is that the popular method due originally to Baumann |jl6[ is adopted: the 
partition function is modified into 

Z = Y.c{N)e-'''N''Z{N) , (5) 

N 

where c{N) = qc^^ for A^ < iVo and c{N) = c^e"^^ for A^ > A'o , while A is 
an adjustable parameter. Of course, one has q/cu = exp (— 2AA'o)- Notice, 
that the parameters c^u need not be specified individually since only their 
ratio enters the detailed balance equations. 

Furthermore, A^ is restricted to A^o =t 6N, with SN <^ Aq- Denote by 
E{N) the experimentally measured density. Since the algorithm satisfies 
the detailed balance equations, within the finite interval under consideration 
E{N) tends to the expression under the sum in (^ for infinite statistics, up 
to global normalization. A sweep of geometry is defined as A"o successive 
steps of the grand-canonical algorithm. 

Another difference, compared to ref. [§, is that we consider Ising spins 
instead of Gaussian fields. When a new vertex is created the new spin(s) 
are generated using the heat bath method. A sweep of geometry is always 



followed by a sweep of spins. The latter is performed using Wolff's |]13 cluster 
algorithm. This algorithm is particularly appropriate in the vicinity of the 
critical spin coupling jSc- Sufficiently above jSc the clusters become large and 
the algorithm loses its efficiency. For such values of /3 the heat bath algorithm 
becomes competitive. We have used it occasionally as an independent check 
of our results. 



Define 



9.„Kiv„) = -i„y^ +« (6) 



Taking the asymptotic form @, one gets: 

E{No + n) ^_2nig-g..)fNo + n 



E{No - n) VA'o - n 



^-2n{g-gcr) "'^ ' (J\ 



and consequently 

9effin,No) = g^r + ^r (8) 

iVo 

to leading order in Nq. The quantity g^ff{n, Nq) is measured for a set of values 
of A^o- Since the right-hand side of (H) is independent of n, one averages over 
n to get an improved estimate of geff^No). The procedure is repeated many 
times and the error on geff{No) is found with the binning method. The 
parameters 7 and gcr are then determined from the linear fit to geff versus 

No'- 

The choice of A is a pure matter of convenience. The experimental his- 
togram E{N) should fall exponentially as one moves off the point N = Nq, 
but one wants to have good statistics for n < 5, say, and therefore this fall 
should not be too fast. The choice of the input value of g involves a subtlety. 
For large A'"o one has roughly 

E{No + n) ^-2n{g-g,r) 



r^u 



-2n(g-gcr) (q\ 

E{No - n) ^^^ 

Thus, the experimental histogram is asymmetric. This asymmetry, together 
with the statistical fluctuations of the histogram, generates a systematic error 
in the measurement of geff{n,NQ). For example, if for a given n one has 
Zexp{No + n) < Zexp{NQ — n), then what is measured is in fact geff{n, Nq) — 5, 
where 6 ~ (l/4n)(e^ — e?_„) and e±„ is the statistical error on E{Nq ± n). 
In practice, one has to start with some exploratory runs in order to get a 
first estimate of gcr, and then use this estimate as the input g in the full 
scale simulation. Proceeding in this way one makes the systematic error 
insignificant compared to the statistical one. 

2.2 Simulating solvable models 

In order to check and gauge our "tools" two solvable models are simulated 
first, viz. pure gravity and the model with 1 spin/site. 

Pure gravity (c = 0): The string susceptibility exponent 7 is known to 
be —1/2. The expression for the partition function is exactly known and the 
critical value of the fugacity z = exp{—g) is Zcr = 27/256 ~ 0.1055. To 
check the code, it is verified that E{N), once normalized, follows very closely 
the exact formula. Then, z is set to value 0.1044, on purpose slightly off 



the critical one, and we perform experiments at A^o=100, 120, 150, 200, 300, 
500 and 1000. The Baumann parameter is given the value A = ln2, forcing 
E{N) to decrease by about a factor 5 between A^o and A^o i 10. The number 
of sweeps per experiment is about 1 to 4x10^ sweeps (following 4000 heating 
sweeps) and the observed acceptance of the grand-canonical algorithm is 
found to be about 72% independently of A^o- 

This exercise is instructive, as it helps appreciating the influence of finite 
size corrections on the results of an experiment of this type. Using all data 
points the estimate 7 = —0.563(33) is produced, significantly different from 
the expected —0.5 . This is not just a fluctuation: the effective exponent 
is below —0.5 precisely because our data follow closely the exact formula. 
However, without the points at A^ < 200 one obtains 7 = —0.509(56), with 
a bigger error (there are less data) but closer to the theoretical value. Of 
course, in this particular model one can determine analytically the threshold 
for the applicability of (^) and the form of correction terms. This is not the 
case for more complicated models. We believe that, while interpreting data 
obtained in simulating such models, guessing the corrections to (H) is a waste 
of time. It is preferable to stick to the simple expression (|) and to accept 
the fact of life that one can only get an effective exponent. Anyway, this is 
the philosophy we have adopted. 

1 spin/site (c = 1/2): The string susceptibility is now 7 = —1/3. The 
critical spin coupling is also known exactly: Per = — | Intanh | In (108/23) ~ 
0.2163 |18|. The code is tested comparing E{N) to the first few terms of 



the expansion in powers of z calculated analytically; the expected values 
are reproduced with an error < 0.1%. We set A = In 2 , 6N = 5 and 
P = Per- The value of fugacity z = 0.05 is found to be acceptable as the 
input one. One sweep is defined as A'o steps of the grand-canonical algorithm 
followed by riy^ = A^/ (cluster size) calls to the Wolff algorithm. The average 
size of clusters produced by the Wolff algorithm is determined prior to each 
production run and it is checked that is remains stable during the run. The 
efficiency of the grand-canonical algorithm is about 57%. It is checked, using 
either a "cold" or a "hot" start, that the system thermalizes after about 2000 
sweeps (cf. ref [jl5|)- This is therefore the number of sweeps to be performed 
initially in order to heat the system. Simulations are carried out for A'q 
= 200, 300, 500 and 1000 with about 1.3 x 10^ sweeps for each A^o- The 
resulting exponent is 7 = —0.317(102). This is unprecise but correct and, 
without further ado, we go over to the truly interesting case of 4 spins/site. 
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2.3 Surfaces with 4 spins/site 

Now, neither the string susceptibihty exponent, nor the critical spin couphng 
Pc are known theoretically. Of course, it is expected that 7(/3) has a maximum 
sX j3 = j3c and falls off to —1/2 as j3 moves away from its critical value. 
The exponent 7 is estimated for several values of (3 in order to observe the 
maximum in question and the corresponding value of 7. We proceed as in the 
1 spin/site case, making the same tests (with similar results) and choosing 
the same values of A and 6N . The fugacity is always chosen so that E{N) is 
almost symmetric in the neighborhood of A^o- The observed efficiency of the 
grand-canonical algorithm is about 56% . The experiment is carried out for 
A/'o=200, 300, 500 and 1000, with 2000 heating sweeps followed by 1.3 x 10^ 
production sweeps, for each A^q- 

The results are shown in fig. 1. The exponent 7 takes its maximum value 
near (3 ^ 0.196. This is reassuring, since this is the critical coupling found 
from microcanonical simulation in ref. (in comparing with their value 
remember that they have been working on the dual lattice). 

The relatively slow fall of 7 on the right of the critical point is associ- 
ated with the worsening performances of the Wolff method (size of clusters 
increasing too much). We have also made some runs with the heat bath 
algorithm. With the latter, in order to get results close to those with Wolff, 
one has (for the system with 4 spins/site) carry out 3 sweeps of spins after 
each sweep of the geometry. And, as expected, the autocorrelation time is 
longer (however not larger than ~ 100 sweeps, this low number is presumably 
due to the fact that the grand-canonical algorithm permanently destroys and 
creates spins). But at /5 = 0.24 the value obtained using the heat bath algo- 
rithm is —0.487(83), which is closer to the expected —1/2 than the rightmost 
point in fig. 1. 

3 DISCUSSION AND SPECULATIONS 

Using the two top values from fig. 1 we get the estimate 7 = —0.195(58). 
It is compatible with the value found in ref. WM, which is —0.167(8). The 



respective errors should not be directly compared, since they have different 
significance. Our error estimate is very conservative. The error given in 



ref. IjTj] is the error obtained in fitting the baby universe distribution in a 



single experiment. Contemplating the dispersion of 7's shown in fig. 1 of ref. 
[Q and corresponding to nearby values of /3, one suspects that this error 
should be multiplied by at least a factor of 2. Nevertheless, our experiment 
is admittedly less precise, simply due to our lower statistics. 

The size of our lattices is the same as the size of baby universes observed 



in ref. |]IJ] . However, with their method one can easier extend simulations to 
larger systems, and it is this which makes it superior to the grand-canonical 
technique employed in our work. Anyhow, it is certainly a good thing to 
have two independent and compatible estimates of the string susceptibihty 
exponent at c = 2, obtained with two completely different methods. 

We wish to end this paper with a speculation concerning c > 1 models. 
We conjecture that the continuum c = 1 theory acts as an "attractor" within 
a finite interval Cq > c > 1. More precisely, for these central charges, we 
suspect the long distance physics of the different models to be the same, 
modulo a renormalization of the length scale. In particular, this means that 
the microcanonical partition function behaves as follows: 

pOcrN 

where A = A(c) is some constant and the dots represent sub leading non- 
universal contributions. One can check that the normalized distribution of 
"baby universes" given in ref. |T^ for c = 2 is compatible with (p!OD, even 
without sub-leading corrections, provided one sets A ~ 20. With such a value 
of A the effective 7 calculated from (|T0|) for N G [200, 1000] is also close to 
the one measured in this work. Such a simple rescaling does not reproduce 
the data |]T^ Q for c > 3. Of course, one could reconcile (|iy) with these 
data playing with the sub-leading terms. This would be a futile exercise, 
however, in view of the lack of any theoretical prediction about the latter. 

Our conjecture is simple and refutable. We obviously expect that the 
effective susceptibility exponent tends to zero when data with larger N are 
used, provided 1 < c < cq. Actually, the authors of ref. []T3[ report that 



for c = 1.6, an extension by the factor of 10 of the size of "baby universes" 
produces a shift of the effective 7 from —0.212(15) to —0.175(24). Because 
of the errors, the effect is not very significant. However, it goes in the right 
direction and the magnitude of the shift, if taken at the face value, is just 
what we expect with A = 10 to 20. This is perhaps not an accident. If true. 



our conjecture also explains why the critical exponents (fractal dimension 
and magnetic exponents) found in microcanonical simulations [^, |l^, 11, 12 
show so little dependence on c just beyond the "barrier". 

The general arguments reviewed in ref. |^ indicate that 7 > implies 
that this exponent has a finite value, which for spin models is most likely | 
(see also ref. [|1^]). We thus expect that, in the continuum limit, the exponent 



7 jumps from the value it takes in the interval [1, Cq) to a moderate value, 
say I , for c > Cq. We imagine that the jump occurs when the density of 
baby universes increases above some critical value. With this perspective, 
we are tempted to consider the small positive values of 7 (e.g. 7 ~ 0.04 to 
0.06 for c = 2) extracted from series expansion in ref. P] as supporting our 
conjecture. 

It has been suggested long time ago (first by David in ref. [^; see also 
[^ and references therein) that the effective 7 becomes positive somewhere 
in the neighborhood of c = 4. The most recent data presented in ref. Wi 



indicate that this happens for 3 ^ c < 4 This result suggests the crude 
estimate 3 < cq < 4. It is perhaps not unreasonable to think that such a 
value of Co is large enough to warrant the approximate validity of the mean 
field calculations for c > Cq. 

Although our conjecture differs from that put forward in ref. ^^ as a 
plausible "c > 1 hypothesis", it is likewise pessimistic: in both scenarios the 
physics of a c > 1 model is reducible to the known theories with c < 1. 

This study has been triggered by a correspondence with D. Johnston, 
whom we wish also to thank for interesting conversations. We wish to thank 
J. Ambj0rn for useful correspondence. One of us (JPK) is indebted to 
the French-British collaboration fund Alliance for covering expenses asso- 
ciated with his visit at the Heriot-Watt University. Laboratoire de Physique 
Theorique et Hautes Energies is Laboratoire associe au C.N.R.S. 
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Figure caption 



Fig. 1 - The string susceptibility exponent 7 versus the spin coupling (3. 
The arrow indicates the critical coupling calculated in ref. [0]. 
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